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The existence of a neutron scattering resonance at a wavevector q* implies a sign change of the gap 
between two Fermi surface regions separated by wavevector q* . For the Fe pnictides, a resonance has 
been observed for a wavevector q* which connects a hole Fermi surface around the F point with an 
electron Fermi surface around the X or Y points of the 1 Fe/unit cell Brillouin zone. Here we study 
the neutron scattering resonance for a five orbital model within an RPA-BCS approximation. Our 
results show that both sign-switched and extended s-wave gaps are consistent with the present data 
for q* near (tt, 0) and that scattering at other momentum transfers can be useful in distinguishing 
between gap structures. 



PACS numbers: 

Recent neutron scattering studiesii^ii^ find a resonance 
in the superconducting state of the 122 Fe pnictides 
which appears below Tc and has a wave vector q* ~ (tt, 0) 
in the unfolded (1 Fe/cell) Brillouin zone. One would 
like to understand what this resonance implies about the 
structure of the superconducting gap A(A;). As is well 
known, the occurrence of resonances in the neutron scat- 
tering spectrum depends through the BCS coherence fac- 
tors on the relative signs of the gap on different parts on 
the Fermi surface separated by q, and thus gives insight 
into the momentum structure of the superconducting 
gap. Nuclear resonance Knight shift measurements^!^ 
support a singlet pairing state and Aig and Big gaps 
have been found in various fluctuation-exchange calcula- 
tionsiiSiiSiiiS as well as renormalization group studiesj^iii^ 
Motivated by the experiments and results from these cal- 
culations, we propose to examine the relationship of the 
resonance to its q* space location and the fc-dependence 
of the gap. Previous calculationsi^ii^ have found that a 
resonance occurs at a wave vector q* ~ (tt, 0) for a sign 
switchedi^ s± gap. This gap has a sign change between 
the inner hole Fermi surfaces that surround the P point 
and the electron Fermi surfaces around the X{tt,0) and 
y(0, tt) points of the unfolded (1 Fe/cell) Brillouin zone. 
One of these calculations was for a two-orbital modeli^ 
and the other used a four-band model but neglected the 
role of orbital-band matrix elements.^'* To adequately de- 
scribe the region of the Fermi surface of the Fe-pnictides 
one needs at least three orbitalsi^ and the orbital-band 
matrix elements are known to play an essential role in 
determining the q dependence of the magnetic suscepti- 
bilityii^ Here we will use an RPA-BCS approximatiortf 
to calculate the neutron scattering response for a five- 
orbital tight binding model with onsite Coulomb and ex- 
change interactions. This model has sufhciently many 
orbitals to describe the electronic structure in the rele- 
vant regions near the Fermi energy. We will also take 



account of the orbital-band matrix elements. 

Figure [1] shows the Fermi surfaces for this 5-orbital 
model with tight binding parameters fit to reproduce the 
LDA bandstructurei^ near the Fermi energy. For a dop- 
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FIG. 1: Fermi surfaces for the 5-orbital model at a doping 
X — 0.125. There are two hole Fermi surfaces ai and a2 
around the F point and two electron Fermi surfaces /3i and 
/32 around the X(7r,0) and y(0,7r) points, respectively. 

ing X = 0.125, there are two hole Fermi surfaces ai and 
a2 around the P point of the unfolded (1 Fe/cell) Bril- 
louin zone and two electron Fermi surfaces /3i and [32 
around the X and Y points. For certain ranges of in- 
teraction parameters, the leading pairing instability in 
an RPA fluctuation exchange approximation occurs for 
an Aig extended s-wave gap and for other regions a Big 
d-wave gap is favored. In the following, we will exam- 
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ine the possibility of a resonance in the inelastic neutron 
scattering for each of these cases as well as for an Aig 
sign-switched s-wave gap. 

The model that we will study consists of a 5-orbital 
(dxz, dyz, dxy, d^2_y2, d^z2_j.2) tight-biuding fit to the 
DFT bandstructure of Cao et alr^ and onsite intra- 
orbital f/, inter-orbital V Coulomb and exchange J inter- 
actions. The Hamiltonian is discussed in the Appendix 
and energies will be measured in units of the largest 
hopping matrix element ty^, i.e. the nearest-neighbor 
hopping along the y-direction between dxz orbitals. For 
this multiorbital problem, the dynamic spin susceptibil- 
ity Xijil:^) that determines the neutron scattering in- 
tensity depends upon an orbital-dependent spin suscep- 
tibility tensor-i2i22, given by 



dTe'--^TrSL^{q,T)Sf{-q,0)) (1) 



Here r and s label orbital indices (1, . . . , 5) corresponding 
to the {dxz, dyz, dxy, d^^-ya, d^z^-r^) orbitals. S'l^iq) = 
5 Sfe d'^-^{k + q)al^fjdsi{k) is the spin flip up operator act- 
ing between the r and s orbitals and S'!I'{—q) = [S'^{q)Y . 
Carrying out the usual analytic continuation of Matsub- 
ara frequencies to the real frequency axis, the dynamic 
spin susceptibility is given by 



E 



Xu{q,iu;m ^ uj + i6) 
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We will approximate the orbital spin susceptibility by 
an RPA-BCS form 



(3) 



with an interaction tensor given in the Appendix. 
This interaction tensor contains onsite intra- and inter- 
orbital Coulomb interaction U and V along with the 
intra-orbital exchange J and pair hopping term J'. 
Here xo is the BCS susceptibility tensor 

(xo)r:(g,t^m) = E M';:^^{k,q) 

X {G"'(fc + g, w„ -f Wm)^ (/c, w„) 

+ F%-k-q,~oj^~oj^)F'''{k,uj,,)lA) 

with 

n^tu \ i^n + E^{k) A(fc) 



Motivated by our previous calculations of the pairing 
interaction and gap functions for the 5-orbital model, we 
will examine the neutron scattering response for a typi- 
cal set of interaction strengths V = U — |J, J — U/8, 
J' = J/2 with [7 = 4 in units of ty^ at a doping x — 0.125. 
The static RPA susceptibility xxl,'^ = 0) the normal 
state is plotted in Fig. [2l Here the (tt, 0) peak of the 
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FIG. 2: The normal state static RPA spin susceptibility 
x[q, 0) normalized to the single particle density of states 
N{Ef) per spin for a doping x = 0.125 with V — U — jJ, 
J — U/8 and U — 4. The peaks are associated with scattering 
between the Qi and Pi and /32 Fermi surfaces shown in Fig.[T] 

undoped system has moved to an incommensurate wave 
vector q* = (tt, O.IStt). This peak is associated with scat- 
tering between the ai and (3i Fermi surfaces. The ridge 
near (tt, tt) is associated with scattering between the /3 
Fermi surfaces. The superconducting gap arising from 
the fluctuation-exchange pairing interaction for this case 
can be parameterized near the Fermi surfaces in terms of 
low order harmonics. For these interaction parameters, 
the leading pairing instability occurs in an Aig channel 
and the next leading instability is in a Big channel. As 
discussed in Ref. [lO|, the pairing strengths for these two 
gap symmetries are quite similar. 

We will consider both of these possibilities using pa- 
rameterized gaps given by 



Ay (A;) — Ay (cos k^ ± cos ky) 



(7) 
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Here the plus sign corresponds to the Aig (extended s- 
wave) gap and the negative sign to the Big (d-wave) gap. 
The amplitudes A^ are adjusted so that the maximum 
magnitude of the gap on the ai and (3 Fermi surfaces is 
0.1 and 0.05 on the a2 Fermi surface. Based on the re- 
sults from the fluctuation-exchange calculation, we have 
^rstui^ 1 q) — '^v (fc + 9)fly'(fc)'it'(^)'^"(^ + • (6) taken Aq2 with the opposite sign to Aq J for the case. 



^ (5) 

and £^{k) = ^jEl{k) + A^{k). The band energies E^{k) 
are measured relative to the Fermi energy and the orbital- 
band matrix elements ay(fc) enter in determining 
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We will also calculate the inelastic neutron scattering for 
the case of a sign-reversed s-wave system in which the 
gaps are isotropic on each Fermi surface with A^^ =0.1, 
Aq,2 = 0.05 and A^^ — Ap^ — —0.1. In each of these 
cases, in order to describe the behavior of the order pa- 
rameter away from the Fermi surface, we have multiplied 
A^(fc) by a Gaussian cutoff exp[-{E{k)/ AE)"^] with AE 
of order several times the gapiS. The resulting gaps are il- 
lustrated in Fig. [3] The anisotropic Aig s-wave gap struc- 
ture in Fig.[3]D is similar to that obtained in both RPAiiiS 
and functional renormalization group calculationsi^. De- 
pending on the interaction and doping parameters, the 
variation of the gap on the (3 Fermi surfaces may or may 
not be sufficient to produce nodes. The cZ-wave Big gap 
illustrated in Fig. [5J: is also similar to the gap structure 
found in these calculations. As noted, for the parameter 
set we have chosen it is the second leading instability but 
with other parameter choices it can become the leading 
instability. 

The inelastic neutron scattering intensity is propor- 
tional to the imaginary part of xil,^)- In Figure 3] we 
show x" ^) versus oj for the normal state and the three 
different superconducting gaps. These calculations have 
all been carried out at a temperature T = 0.005 and 
the normal state has simply been included for compari- 
son purposes. Figure shows results for a momentum 
transfer qi ~ (tt, 0.157r), where we expect from Fig. [T] 
that a resonance will appear. When the gap opens in 
the superconducting state, low energy spectral weight 
is shifted to higher energies and the quasi-particle-hole 
damping is suppressed. The resonance appears when the 
real part of the denominator of Eq. [3] vanishes provided 
that ReA{k)A{k + q*) < so that the coherence factor 

If A{k)A{k + q*) \ 

2 V £{k)£{k + q*) J ^ ' 

goes to 1 rather than zero. Thus, just as in the well 
known cuprate case, the observation of a resonance pro- 
vides evidence of a change in the relative signs of the 
gap between two regions of the Fermi surface^^'^^. Here 
for qi = (tt, 0.157r) shown in Fig. [3^, the resonance im- 
plies that this sign change occurs between regions on the 
ai and (5 Fermi surfaces. For the sign switched s± gap 
shown in Fig. [3^ this sign change occurs over the entire 
Fermi suface. For the extended s-wave gap illustrated 
in Fig. the qi wavevector again dominantly connects 
regions of the a and (3 Fermi sheets which have opposite 
signs. Because of the gap nodes and gap anisotropy on 
the (3 Fermi surface sheet, the response is weaker than 
that for the sign-switched s-wave. 

For q2 = (tt, 0.567r) shown in Fig. [3)d, there are im- 
portant contributions involving particle-hole scattering 
between the (3i and /?2 Fermi surfaces. In this case, the 
sign switched s-wave gap has the same sign on both of 
the P Fermi surfaces so that the coherence factor, Eq. [8] 
vanishes. Then as shown in Fig. Uja, there is no reso- 
nance for the gap and one only sees that spectral 
weight is shifted from low frequencies in the normal state 
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FIG. 3: Plots of the gap A^(fc) for the ai, a2, and (3 Fermi 
surfaces (a) the sign-switched s-wave, (b) the A\g (extended 
s-wave) gap Ay(cosfca; -|-cosfcy) and (c) the Big (d-wave) gap 
Ai, (cos kx — cos ky). All are plotted with the Gaussian cut-off. 
Here red denotes positive and blue negative values of the gap. 
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FIG. 4: The RPA-BCS dynamic spin susceptibility % il*,^) 
versus w for (a) q* = (tt, O.IStt), (b) q* = (7r,0.567r) and (c) 
q = (0.727r, 0.727r) normalized to N{Ef). The two Aig gaps 
exhibit a resonance and the d^2_y2 gap has a broad response 
at (7r,0.157r). Only the Big d-wave and the Aig extended s- 
wave show a response at (tt, O.SGtt) and at (0.727r, 0.727r) only 
the d^2_y2 gap has a resonance. 



to frequencies above ^ 2A in the superconducting state. 
However, the resonance is present for both the extended 
s-wave state and the Big (d-wave) state since, as seen 
in Fig. [21 92 connects regions of the f3i and /32 Fermi 
surfaces where there is a sign change of the gap and 
ReA{k)A{k + q*) < 0. Again these regions are smaller 
for the extended s-wave. Fig. [3)3, compared to the dx2_y2 
gap shown in Fig. [3j:. 

For the ^3 = (0.727r, 0.727r) momentum transfer shown 
in Fig. [Sjj, only the d^^-y^ gap leads to a resonance re- 
sponse. In this case, as in the previous case, the scat- 
tering involves transitions between the /3 sheets. For 
the da;2_y2-wave gap. Fig. [31:, these sheets have opposite 
signs, but for the extended s-wave, Fig.[3}3, these regions 
have the same sign and the resonance is suppressed. 

These calculations show that both the sign-switched 
s-wave and the extended s-wave gaps are consistent with 
the occurrence of a resonance in the neutron scatter- 
ing response which is observed in the superconducting 
state of the 122 Fe pnictides near (tt, 0) in the un- 
folded (1-Fe/zone) Brillouin zone. In our calculations, 
this resonance is strongest for the phenomological, sign- 
switched gap. However it also appears as a clear reso- 
nance for an extended s-wave gap and even as a weak 
broad structure for a dx2_y2 gap. At a momentum trans- 
fer 92 = (tt, O.SGtt), the transitions involve scattering be- 
tween regions of the (ii and /32 Fermi surfaces where the 
sign-switched s-wave gap has the same sign. For the di- 
agonal momentum transfer shown in Fig. [H;, only the 
d-wave gap exhibits a resonance. Thus further neutron 
measurements at other momentum transfers can narrow 
the possible gap structure. Our results suggest that cuts 
in q-space going from (tt, 0) to (tt, tt) and along the diag- 
onal direction in the unfolded 1-Fe zone are most promis- 
ing in providing useful information. 
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I. APPENDIX 

The full DFT band structure can be fitted in the vicin- 
ity of the Fermi energy using a tight-binding approxima- 
tion with the five Fe d orbitals as basis set. Here we intro- 
duce a coordinate system aligned parallel to the nearest 
neighbor Fe-Fe direction. The Hamiltonian for this 5 or- 
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bital model can be written as 

<^rnSmn)dl^^{k)dna{k) (A-1) 

k,a mn 

where d\^^{k) creates a particle with momentum k and 
spin (T in orbital m. The symmetry of the kinetic energy 
terms {k) can be derived from a Slater-Koster based 
parametrization that respects the symmetry of the FeAs 
layers. The exact sizes of the hopping parameters can be 
used as fitting parameters to approximately reproduce 
the band energies at the high symmetry points (F, X, 



M) and also the band structure along the high symmetry 
directions - at least in the vicinity of the Fermi energy. 
The explicit form of the kinetic energy terms S^mn as well 
as the on-site energies e„i and the values of the hopping 
terms used for the fitting of the band structure by Cao 
et alM- can be found in the appendix of Rcf. 9. 
The interaction tensor U" is given by 

(UIZ = u , mtS = J/2 , (uTat = ^+vAun^l - J' 

(A-2) 

where a ^ b. 



Electronic address 
Electronic address 
Electronic address 
Electronic address 



: "maierta@ornl.gov' 
: graser@phys.ufl.e duj 
: djs@physics.ucs b.edu| 
: pjh@phys.ufl.edu' 
A.D. Christiansen et al, arXiv:0807.3932 (2008). 
M.D. Lumsden et al, arXiv:0811.4755 
S. Li et al, preprint arXiv:0902.0813 (2008). 
H.-J. Grafe et al, Phys. Rev. Lett. 101, 47003 (2008). 
K. Ahilan et al, arXiv:0804.4026 
Ning et al, arXiv:0808.1420 

Kurok i et al, Phys. Re v. Lett. 101, 87004 (2008). 
Ikeda, larXiv:0810.1828 i 

Yanagi, Y. Yamakawa and Y. One, ' arXiv:0808."TT92] 
0809.3189 

S. Graser, T.A. Maier, P.J. Hirschfeld and D.J. Scalapino, 
New J. Phys. 11, 025016 (2009). 

A.V. Chubukov, D.V. Efremov and I. Eremin, Phys. Rev. 



B 78, 134512 (2008). 

^2 F. Wang et al, arXiv:0807.0498] 

T.A. Maier and D.J. Scalapino, Phys. Rev. B 78, 20514(R) 
(2008). 

M.M. Korshunov and I. Eremin, EPL 83, 67003 (2008). 
I.L. Mazin et al, Phys. Rev. Lett. 101, 57003 (2008). 
P. A. Lee and X.-G. Wen, Phys. Rev. B 78, 144517 (2008). 
" N. Bulut and D.J. Scalapino, Phys. Rev. B 47, 3419 (1993). 
C. Cao, P.J. Hirschfeld and H.-P. Cheng, Phys. Rev. B 77, 
220506(R) (2008). 

T. Takimoteo, T. Hotta and K. Ueda, Phys. Rev. B 69, 

104504 (2004). 

^° T. Kariyado and M. Ogata, preprint arXiv:0812.4664l 
(2008). 

N. Bulut and D.J. Scalapino, Phys. Rev. B 53, 5149 (1996). 
^2 M. Eschrig, Adv. Phys. 55, 47-183 (2006). 



